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PART I This part carries 40% of the marks for the examination. we 


Answer ALL the eight questions in this part. Not all the questions in this part carry. 
egual marks. 


Question 1 If F(x, y) is a function defined on the (x, y)-plane, and (u, v) is a coordinate system 
in this plane defined by $ 4 


/ 


S 


x=3u? +0? and y= —w, 


OF. OF (OF 

express zz in terms of u, v, Ju and a: 

Question 2 Laplace’s equation in two-dimensional polar coordinates (r, 6) is 

Ou lôu, 1 
+ art ao 

The method of separation of variables yields as a solution 
u = (Ar* + Br-*)(Ccosk@ + D sin k0) , 

where A, B, C, D and k are constants, and k > 0. 


=0 (r>0,all6). 


Write down: 
(i) the most general solution of this type for 
0<r<a, all@ 
satisfying 
u(a, 0) = 0 (all 6) 
u(r, 0) = u(r, 0 + 2n) (<r«<a), 
where a is a positive constant; 
(ii) the most general solution for r > a, 0 < 0 < n satisfying 
lim u(r, 0) = 0 (0<6<n) 
=o 


Aners (r >a). 


Question3 If S is a surface in 3-dimensional space bounding a region R, and v is a function 
dv 


ôn 


use the Divergence Theorem to prove that v = 0 in R and on S. 


that is harmonic in R and satisfies — + kv = 0 on S, where k is a positive constant, 


You may use the result 


div (u grad v) = uV?v + (grad u) - (grad v) . 


Question 4 The explicit scheme 
1 2 
k Atij = R luig + yy 
is to be used to obtain an approximate solution of the parabolic equation 
U(x, t) _ 2 U t) 
a Ox? 


Show that the order of the local truncation error is O(h?) + O(k), where h and k 
are the step sizes used for x and ¢ respectively. 


+ U(x,t). $ 
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Question 5 


Question 6 


Question 7 


Either 7(a) 


Or 1 (b) 


Consider the problem 
aux 
ME) _ (qcos x ~ aux) = 0 (0<x<3) 
u(0) = u (5) =0, 


where q is a positive constant and 2 is an eigenvalue. Using a comparison theorem, 
show that the lowest eigenvalue lies between 4 and 4 + q. 


(i) Which of the following equations has two distinct families of characteristics? 


au eu 

at Se FF =P 0) 
ĝu ĝu 2 ou 

aa a ae = 9 @) 


(ii) Consider a problem with initial data on the line ¢ = 0 for the equation you 
chose in part (i). Draw a sketch showing the domain of dependence of the 
point where x = 1, £ = 3, giving the equations of any lines or curves you draw. 


Answer either 7(a) or 7(6). Do NOT attempt both. 


If u(x) satisfies the differential equation 
d?u(x) 
dx? 
show by putting £ = e~* that the corresponding function u(t) satisfies 


2, 
ÈO AO, naa = 0 O<t<}). (2) 


e 


+u(x)=0 (x>0), (1) 


Hence find the solution of equation (1), given that u(x) has the limit 2 for large x. 
A finite-difference scheme for the equation 


GU(x, t) 3 U(x, t) 


2 
oe UG) ĝt? 


is 

Uige T Yi F r(t; j+)? Gisiger — Mises + Una yjei) 0 
where x = ih, t = jk, and r = k/h’. It is required to use this scheme to solve the 
equation with the boundary condition U(x, t) = 0 when x = 0 and when x = 1, 
using a mesh with h = $. (See diagram.) 


The values of u at the points A, B are supposed known (call them u4, ug). 
(i) Write down the pair of equations which must be solved to find uc and up. 


Gii) Show how to solve these equations by Newton’s iterative method, and comment 
on its convergence. 


Turn over 
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Question 8 Briefly describe in your own words the concepts of consistency, convergence, and 
stability in connection with finite-difference approximations to differential equations. 


PART 0 


Question 9 


Question 10 


Question 11 


This part carries 60%, of the marks for the examination. 
Answer at most FOUR questions from this part. All questions carry equal marks. 


(i) In which parts of the region {(x, 1):¢ > 0} in the (x, ¢)-plane is the equation 
oy Fu 
oxi  * g 


elliptic, in which parts is it hyperbolic, and in which parts is it parabolic? 


1 =0 (@>0 


(ii) In the region (call it R) where this equation is hyperbolic, find the 
characteristics. 


(ii) Briefly explain the importance of characteristics. 
(iv) Find a coordinate system in R for which the equation takes its standard form. 


(v) Carry out the transformation to standard form. 


An inflatable tent is made by attaching the edges of a rectangular sheet of suitable 
material to the ground and then pumping in air underneath. The height of the 
material above the ground at the point (x, y), which we denote by u(x, y), satisfies 
the equation 


Vu=—-F, 0) 


where F is a positive constant depending on the excess air pressure inside the tent. 
If the lengths of the sides of the rectangle are a and b, write down the boundary 
conditions for the equation (1) in a coordinate system whose axes are along two 
sides of the rectangle. 


Write down a formula giving the volume of air inside the tent in terms of the function 
u(x, y). 

Using an extremum principle, write down an expression which gives a positive 
lower bound on this volume. Evaluate this lower bound for the case when a = b. 


Describe one explicit and one implicit finite-difference scheme for numerically 
solving the problem 


, . 
Wen 3G) O<x<1,t>0) 
UO, t) = 0 

boro 
UU, 0) =0 


Ux, =f» Ox<x<1) 
where f(x) is a given function. 


Compare and contrast your two schemes with respect to stability, accuracy, and 
ease of computation, giving reasons (but not detailed proofs) for any assertions 
you make. 


Turn over 
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Question 12 


Question 13 


@ 


Gi) 


Consider the problem 


U(x, t) 4 U(x, t) 
a ax? 


ôU, 1) _ o) 
ôx fort >0 
Ud, t) = 0 
U(x, 0) = f(x) Uers; 


where f must satisfy suitable conditions but is otherwise arbitrary. Using the 
method of separation of variables, find the formal Fourier series representing 
U(x, t) in terms of f. 


(0<x<1,t>0) 


Consider next the Fourier series found by you in part (i), evaluated at ¢ = 0. 
Suppose the functions (a) and (b) below are each substituted for the initial 
function f in this series. 


For each function state at which points x, with 0 < x < 1, the resulting series 
does not converge pointwise to U(x, 0). If there are such points, state the value 
of the sum of the series there. 


Note carefully that no integrations are necessary in order to answer this question. 


s=) 
@ fx) =1-x O<x<1) s 
fx) 
x @<x<}) 
orre= a Glen x 


(iii) The function (a) of (ii) above satisfies f'(x) = —1 for all x in its domain, 


(i) 


ii) 


yet if the Fourier series for this function, as used in (i) above, is differentiated 
term by term, the resulting series sums to zero when x = 0. Explain briefly 
how the concept of uniform convergence helps to resolve this apparent 
contradiction. 


Write in self-adjoint form the equation 


d d 
O 2g MO) + yx) = FO), (1) 


where æ is an unspecified number and F an unspecified real-valued function. 


Show that if |aJ < 1 the influence function of the self-adjoint equation is 
given by 


R(x, č) = w~! e**® sin w(x — €) , 


where w = /(1 — a”). 


(iii) Obtain the solution of (1) which satisfies the conditions 


u0) =0, vO) =0, 


expressing your solution in terms of R. 


(iv) Show, if you have not already done so earlier in the question, that R(x, ¢) and 


(2 


its derivatives have the properties characterizing the influence function, and 
that the solution obtained in (iii) is linear in F. 


Give a physical interpretation of R(x, ¢). 
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Question 14 (i) Find a 5-point formula for the solution of 


UC, U(x, 
E e mo o 


on a square mesh of side h. 


(ii) It is desired to solve equation (1) on the rectangle 
0<x<l, O<y<3, 
with the boundary conditions 
UO,» = U, y)=0 O<y<% 


U(x, #) = 0 
Uo- $2 (x, 0) = poers 1). 


Write down the system of equations for the approximate values of U at the 
four points marked A, B, C, D in the diagram, arising when the 5-point 
formula is used, with central differences for dU/dy, and h = }. 


(iii) Write down the iteration matrix for Jacobi’s method of solving this system of 
equations and state whether this method is convergent, giving the reason for 
your answer. 


(iv) Is the S.O.R. method convergent for this system of equations? Give your 
reason. 


Turn over 
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Question 15 


Question 16 


Consider the Sturm-Liouville problem 


5 (oe E — g(x) ux) + A p(x) ux) =0 (0<x<1) 


u(0) = (1) = 


where p and q are continuous, p is continuously differentiable, p and p are positive, 
and q is non-negative, for 0 < x < 1. 


_ (i) Show that, if u(x) is a non-trivial solution of the integral equation 


1 
ula) = Af Gx, DDD, 
where G is the Green’s function for the problem 


ál p(x) 2a) - g(x) ux) = - F(x) 


u(0) = ul) = 0, 


then u(x) is an eigenfunction, and A is an eigenvalue, of the Sturm-Liouville 
problem. Show all steps in your argument and state what theorems you are 
using. 


(ii) Prove, stating what theorems you are using, that 


f f, Gtx D ols) E de dk = > T 


where A,, Àz» . . . are the eigenvalues of the Sturm-Liouville problem. 


A non-uniform string is stretched with tension T between the points x = 1, x = 2 
on the x-axis of a rectangular coordinate system, both ends being held fixed. 
The density of the string at the point x is m/x?, where m is a positive constant. 


(i) Find the differential equation and boundary conditions that model small 
transverse vibrations of the string, if gravity is ignored. 


(ii) Show that your differential equation has solutions of the form 
el! f(x) , 
where f satisfies 


z8? 
2," mar fis) = 0. 


(iii) By means of the substitution x = e*, or otherwise, show that the angular 
frequencies, w, are given by 


amo? _ | 2nn\? i3 
T = + \in2 (n= 1,2,...). 
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